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Abstract — This  paper  proposes  a  new  analysis  method  of 
multi-state  systems  with  multi-state  components  using  multi¬ 
valued  decision  diagrams  (MDDs).  The  multi-state  systems  with 
multi-state  components  can  be  considered  as  multi-valued  func¬ 
tions,  called  structure  functions.  Since  the  structure  functions 
are  usually  monotone  increasing  functions,  they  can  be  rep¬ 
resented  compactly  using  edge-valued  MDDs  (EVMDDs).  This 
paper  proposes  an  efficient  analysis  method  using  EVMDDs. 
It  shows  that  by  using  EVMDDs,  the  structure  functions  can 
be  represented  more  compactly  than  existing  methods  using 
ordinary  MDDs,  and  systems  can  be  analyzed  with  comparable 
computation  time. 
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I.  Introduction 

In  recent  years,  various  systems  such  as  computer  server 
systems,  telecommunication  systems,  water,  gas,  and  elec¬ 
tricity  distribution  systems,  become  tolerant  to  faults  and 
errors.  Even  if  a  fault  occurs,  these  systems  still  keep  on 
working  with  an  acceptable  or  degraded  performance  level. 
Thus,  these  systems  cannot  be  modeled  by  two  states:  work¬ 
ing  and  failure.  In  addition,  with  the  advance  of  technology, 
each  component  in  a  system  also  becomes  fault  tolerant.  To 
model  such  systems,  a  multi-state  system  with  multi-state 
components  are  often  used  [11],  [14],  [16]. 

Fault  tolerance  is  usually  achieved  by  multiplexing  com¬ 
ponents  in  a  system.  However,  multiplexing  noncritical 
components  to  design  a  fault  tolerant  system  is  inefficient 
and  not  cost  effective.  Also,  if  critical  components  are  not 
sufficiently  tolerant  to  faults,  fault  tolerance  of  the  system 
is  not  sufficient.  Thus,  identifying  which  components  are 
critical  to  achieve  fault  tolerance  of  the  system,  and  mul¬ 
tiplexing  them  are  important,  especially  for  safety-critical 
systems  such  as  flight  control  systems  and  nuclear  power 
plant  monitoring  systems  [2], 

To  identify  critical  components  and  system  weaknesses, 
analyzing  multi-state  systems  by  various  assessment  mea¬ 
sures  is  required  [11].  Among  them,  assessing  the  proba¬ 
bility  of  each  state  of  a  multi-state  system  is  essential  to 
the  design  of  a  dependable  fault  tolerant  system  [14],  [16]. 
Various  methods  to  analyze  multi-state  systems  efficiently 
have  been  proposed.  Many  existing  methods  are  based  on 
the  Markov  model  [3].  However,  they  are  impractical  for 


a  large  multi-state  system,  since  their  time  complexity  is 
0(m?n),  where  m  is  the  number  of  states,  and  n  is  the 
number  of  components  in  a  multi-state  system  [2].  To 
analyze  large  multi-state  systems  efficiently,  methods  based 
on  binary  decision  diagrams  (BDDs)  [1],  [2],  [4],  [16]  and 
multi-valued  decision  diagrams  (MDDs)  [6],  [13],  [14]  have 
attracted  much  attention. 

Since  multi-state  systems  with  multi-state  components  can 
be  considered  as  multi-valued  functions,  called  structure 
functions,  they  can  be  represented  by  BDDs  and  MDDs. 
Probabilities  of  states  can  be  computed  using  BDDs  and 
MDDs  with  the  time  complexity  proportional  to  the  number 
of  nodes  in  a  decision  diagram.  BDDs  represent  structure 
functions  by  converting  multi-valued  variables  and  function 
values  into  binary  vectors  using  one-hot  encoding  [16].  By 
using  BDDs,  various  analysis  methods  well-established  for 
binary-state  systems  can  be  directly  applied  to  multi-state 
systems.  However,  converting  into  binary  vectors  produces 
many  binary  variables  and  many  binary  functions,  and 
results  in  large  BDDs.  Thus,  use  of  MDDs  is  more  natural 
and  more  promising  for  larger  multi-state  systems. 

For  recent  large  and  complex  multi-state  systems,  how¬ 
ever,  decision  diagrams  that  represent  systems  more  com¬ 
pactly  are  desired.  Since  structure  functions  are  usually 
monotone  increasing  [13],  they  can  be  represented  com¬ 
pactly  using  edge-valued  MDDs  (EVMDDs)  [10].  How¬ 
ever,  analysis  of  multi-state  systems  using  EVMDDs  is  not 
straightforward.  As  far  as  we  know,  an  analysis  method 
using  EVMDDs  has  never  been  reported.  Thus,  in  this  paper, 
we  propose  an  efficient  analysis  method  using  EVMDDs. 

This  paper  is  organized  as  follows:  Section  II  defines 
multi-state  systems,  and  EVMDDs.  Section  III  shows  repre¬ 
sentations  of  structure  functions  using  MDDs  and  EVMDDs, 
and  in  Section  IV,  we  propose  an  analysis  method  using 
EVMDDs.  Experimental  results  are  shown  in  Section  V. 

II.  Preliminaries 

This  section  defines  multi-state  systems,  structure  func¬ 
tions,  and  MDDs  to  represent  structure  functions. 

A.  Multi-State  Systems  and  Structure  Functions 

Definition  1:  A  multi-state  system  is  a  model  of  systems, 
which  represents  performance,  capacity,  or  reliability  levels 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

MAY  2012 


2.  REPORT  TYPE 


4.  TITLE  AND  SUBTITLE 

Analysis  of  Multi-State  Systems  with  Multi-State  Components  Using 
EYMDDs 

6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School, Dept,  of  Electr.  and  Comp. 
Eng., Monterey, CA, 93943 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


3.  DATES  COVERED 

00-00-2012  to  00-00-2012 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

Proceedings  of  the  42nd  International  Symposium  on  Multiple-Valued  Logic,  Victoria,  Canada,  May  14-16, 
2012,  pp.  122-127 

14.  ABSTRACT 

This  paper  proposes  a  new  analysis  method  of  multi-state  systems  with  multi-state  components  using 
multivalued  decision  diagrams  (MDDs).  The  multi-state  systems  with  multi-state  components  can  be 
considered  as  multi-valued  functions  called  structure  functions.  Since  the  structure  functions  are  usually 
monotone  increasing  functions,  they  can  be  represented  compactly  using  edge-valued  MDDs  (EVMDDs). 
This  paper  proposes  an  efficient  analysis  method  using  EVMDDs.  It  shows  that  by  using  EVMDDs,  the 
structure  functions  can  be  represented  more  compactly  than  existing  methods  using  ordinary  MDDs,  and 
systems  can  be  analyzed  with  comparable  computation  time. 

15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF 

ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

6 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


X\  X2  X3  X4  f 

~0  0  0  0  0 

0  0  0  1  0 


(a)  Multi-state  system.  (b)  Structure  function. 


Figure  1.  Multi-state  system  for  network  flow  and  its  structure  function. 


of  the  systems  as  states.  It  usually  has  more  than  two  states. 
When  components  in  a  system  are  modeled  as  well,  it  is 
called  a  multi-state  system  with  multi-state  components. 
In  this  paper,  it  is  simply  called  multi-state  system. 

Definition  2:  A  state  of  a  multi-state  system  depends  only 
on  states  of  components  in  the  system.  The  system  with  n 
components  can  be  considered  as  a  multi-valued  function 
f(x i,X2,...,xn)  :  R\  x  R2  x  . . .  x  R„  —> >  M,  where  each  x, 
represents  a  component  with  r,  states,  R,  =  {0, 1 , . . . ,  r,-  —  1} 
is  a  set  of  the  states,  and  M  =  {0, 1, ...  ,m  —  1}  is  a  set  of 
the  m  system  states.  This  multi-valued  function  is  called  a 
structure  function  of  the  multi-state  system. 

Definition  3:  A  multi-valued  function  fix i,X2,  ■  ■  ■  ,x„)  is 
a  monotone  increasing  function  iff  for  any  Xj, 

f(x I,x2,  ■  ■  ■ , Xj~  1 , a,xi+i , . . .  ,xn ) 

<  f(x  1  ,X2, . . .  ,X,'—  1 ,  P,X;+l , . . .  ,Xn), 

where  a,  p  G  R„  and  a  <  p. 

In  many  applications,  states  of  a  system  and  its  compo¬ 
nents  are  totally  ordered,  and  deterioration  of  a  component  in 
the  system  affects  deterioration  of  the  whole  system.  Thus, 
structure  functions  usually  become  monotone  increasing 
functions  by  assigning  a  value  to  each  state  in  ascending 
order  (i.e.  the  worst  state  is  0  and  the  best  state  is  m  —  1  or 

n  - 1). 

Example  1:  Fig.  1(a)  shows  a  multi-state  system  for  net¬ 
work  flow  such  as  water,  gas,  and  electricity  distribution 
systems.  In  this  figure,  each  edge  xt  has  three  states  which 
correspond  to  transmission  capacities:  0  unit  (disconnected), 
3  units  (deteriorated),  and  5  units  (fully  transmittable).  And, 
the  system  has  six  states  which  correspond  to  the  maximum 
number  of  units  that  the  target  node  T  can  receive  from  the 
source  node  S:  0,  3,  5,  6,  8,  and  10. 

By  assigning  six  values  (0,  1,  2,  3,  4,  and  5)  to  these 
states  in  ascending  order,  we  obtain  the  6-valued  structure 
function  /  shown  in  Fig.  1(b).  This  is  a  monotone  increasing 
function.  (End  of  Example) 


Figure  2.  MDD  for  the  structure  function. 


Figure  3.  EVMDD  for  the  structure  function. 


II.  Multi-Valued  Decision  Diagrams 

Definition  4:  A  multi-valued  decision  diagram  (MDD) 

is  a  rooted  DAG  representing  a  multi-valued  function  /. 
The  MDD  is  obtained  by  repeatedly  applying  the  Shannon 
expansion  to  the  multi-valued  function  [5].  It  consists  of  non¬ 
terminal  nodes  representing  sub-functions  obtained  from  / 
by  assigning  values  to  certain  variables.  It  also  has  terminal 
nodes  representing  function  values.  Each  non-terminal  node 
has  multiple  outgoing  edges  that  correspond  to  the  values  of 
multi-valued  variable.  The  MDD  is  ordered;  i.e.,  the  order 
of  variables  along  any  path  from  the  root  node  to  a  terminal 
node  is  the  same.  When  an  MDD  represents  a  function  for 
which  multi-valued  variables  have  different  domains,  it  is  a 
heterogeneous  MDD  [8].  In  the  following,  a  heterogeneous 
MDD  is  also  denoted  by  the  MDD. 

Definition  5:  An  edge-valued  MDD  (EYMDD)  [10]  is 
an  extension  of  the  MDD,  and  represents  a  multi-valued 
function.  It  consists  of  one  terminal  node  representing  0  and 
non-terminal  nodes  with  edges  having  integer  weights;  fl¬ 
edges  always  have  zero  weights.  In  an  EVMDD,  the  function 
value  is  represented  as  a  sum  of  weights  for  edges  traversed 
from  the  root  node  to  the  terminal  node. 

Example  2:  Fig.  2  and  Fig.  3  show  an  ordinary  MDD 


and  an  EVMDD  for  the  structure  function  of  Example  1. 

(End  of  Example) 

III.  MDDS  AND  EVMDDS  FOR  STRUCTURE  FUNCTIONS 

This  section  derives  upper  bounds  on  the  number  of  nodes 
in  an  MDD  and  an  EVMDD  for  a  structure  function.  For 
simplicity,  in  the  following  theorems,  we  assume  that  all 
components  x,  in  a  system  have  the  same  number  r  of 
states  (i.e.,  all  variables  x,  have  the  same  domain).  However, 
generalization  to  a  case  where  all  variables  x,  have  different 
domains  is  straightforward. 

Theorem  1:  For  a  structure  function,  the  number  of  nodes 
in  an  MDD  is  at  most 

i”~l  -  1  / 

- 7~+m  , 

r—  1 

where  /  is  the  largest  integer  satisfying  r"_/  >  mr> ,  m  is  the 
number  of  system  states,  n  is  the  number  of  components, 
and  r  is  the  number  of  component  states. 

(Proof)  See  Appendix. 

Theorem  1  shows  that  the  upper  bound  for  an  MDD 
depends  only  on  m,  n,  and  r.  It  is  independent  of  mono¬ 
tonicity  of  structure  functions.  However,  in  many  applica¬ 
tions,  structure  functions  are  usually  monotone  increasing. 
Thus,  decision  diagrams  suitable  for  monotone  functions 
are  preferable.  Since  EVMDDs  can  represent  monotone 
functions  compactly,  EVMDDs  are  preferable  for  many 
monotone  structure  functions. 

Definition  6:  Let  No  be  the  set  of  nonnegative  integers, 
and  let  p  be  an  integer.  An  integer  function  f(X)  :  No  — >  Z 
such  that  0  <  f(X  +  1 )  —  f(X )  <  p  and  /( 0)  =  0  is  an  M/> 
monotone  increasing  function  on  No,  where  Z  is  the  set  of 
integers.  That  is,  an  M/;-monotone  increasing  function  f(X ) 
satisfies  /( 0)  =  0,  and  the  increment  of  X  by  one  increases 
the  value  of  f(X)  by  at  most  p. 

A  monotone  multi-valued  function  can  be  converted  into 
an  Mp-monotone  increasing  function  by  considering  the  set 
of  multi-valued  variables  x,  as  an  /‘-valued  vector: 

X  =  (xniXn—l , . . .  ,Xl  )>-, 

and  EVMDDs  for  monotone  multi-valued  functions  have  the 
same  complexity  as  EVMDDs  for  M/j-monotone  increasing 
functions  [10].  Thus,  we  derive  an  upper  bound  of  an 
EVMDD  for  an  M //-monotone  increasing  function  in  the 
following: 

Theorem  2:  For  an  M/>monotone  increasing  function 
converted  from  a  multi-valued  function,  the  number  of  nodes 
in  an  EVMDD  is  at  most 

_  i  / 

j  +S(P+  1)'J_1 

/—  1 

where  l  is  the  largest  integer  satisfying  P1-1  >  ( p+  l)^-1, 
and  n  is  the  number  of  variables  in  the  original  multi-valued 
function. 


(Proof)  This  is  the  straightforward  generalization  of  the 
theorem  for  edge-valued  binary  decision  diagrams  (EVB- 
DDs)  [9]  (i.e.,  r  =  2).  Thus,  we  can  also  extend  the  proof 
for  EVBDDs  to  EVMDDs  trivially.  I 

Theorem  2  shows  that  the  upper  bound  for  an  EVMDD 
depends  on  the  value  of  p,  not  on  the  number  of  system 
states  m.  Thus,  even  if  m  is  large,  EVMDDs  have  a  small 
number  of  nodes  when  the  value  of  p  is  small.  In  the  future, 
systems  will  become  more  complex,  and  thus,  m  will  become 
larger.  For  such  systems,  MDDs  require  many  nodes.  On 
the  other  hand,  EVMDDs  can  represent  even  such  systems 
compactly  if  the  value  of  p  is  small. 

IV.  Analysis  Methods  Using  MDDs  and  EVMDDs 

This  section  formulates  a  problem  of  system  analysis, 
and  then  proposes  an  algorithm  to  solve  the  problem  using 
EVMDDs. 

Definition  7:  The  probability  that  a  structure  function 
/  has  the  value  s  is  denoted  by  Ps(f  =  s),  where  s  £ 

{0,1 _ ,m—  1 } .  The  probability  that  a  component  x,  has  the 

value  c  is  denoted  by  Pc  (x,  =  c) ,  where  c  £  { 0 , 1 , . . . ,  n  —  1 } . 

Problem  1:  Given  a  structure  function  /  of  a  multi- state 
system  and  the  probability  of  each  state  of  each  component 
in  the  system  Pc(xj  =  c),  compute  the  probability  of  each 
state  of  the  multi-state  system  Ps(f  =  s). 

In  this  problem,  we  assume  that  all  components  are 
independent  of  each  other.  That  is,  each  state  of  a  component 
appears  independently  of  the  states  of  other  components. 

A.  Analysis  Method  Using  MDDs 

Problem  1  can  be  solved  efficiently  using  node  traversing 
probabilities  in  an  MDD  that  are  introduced  to  compute  the 
average  path  length  on  an  MDD  [7], 

Definition  8:  In  an  MDD,  a  sequence  of  edges  and  nodes 
leading  from  the  root  node  to  a  terminal  node  is  a  path.  The 
node  traversing  probability,  denoted  by  NT P(vj),  is  the 
probability  that  an  assignment  of  values  to  variables  selects 
a  path  that  includes  the  node  v,. 

Since  terminal  nodes  of  an  MDD  for  a  structure  function 
represent  system  states,  node  traversing  probabilities  of 
terminal  nodes  correspond  to  the  probabilities  of  system 
states.  The  node  traversing  probabilities  can  be  computed  by 
visiting  each  node  only  once  in  the  breadth  first  order  from 
the  root  node.  Thus,  the  time  complexity  of  this  analysis 
method  is  0(Nm),  where  Nm  is  the  number  of  nodes  in 
an  MDD.  Other  existing  methods  whose  time  complexity  is 
0(Nm)  also  analyze  multi-state  systems  in  a  similar  way  [6], 
[13],  [14], 

Example  3:  Let  us  compute  node  traversing  probabil¬ 
ities  for  the  MDD  in  Fig.  2.  In  this  example,  we  as¬ 
sume  that  all  states  of  each  component  occur  with  the 
same  probability  1/3.  First,  we  have  NTP(v\ )  =  1  for 
the  root  node  vi  since  the  root  node  occurs  in  all 
paths.  Then,  we  compute  NTP(yi)  =  NTP[yi)  x  1/3  and 


Input:  An  EVMDD  for  a  structure  function  of  a  multi-state  system,  and  the  probability  of  each  state  of  each 

component  in  the  system  Pc(xj  =  c). 

Output:  Probability  of  each  state  of  the  multi-state  system  Ps(f  =  s). 

Step:  The  following  procedures  are  applied  to  each  node  recursively  from  the  root  node. 

1 .  If  a  node  v  has  been  already  visited,  then  return  probabilities  for  the  sub-function  f,  that  have  been  already 
computed.  Else,  go  to  the  next  step. 

2.  If  the  node  v  is  the  terminal  node  T,  then  return  the  probability  for  the  constant  zero  function:  Ps(fr  =  0)  =  1. 
Else,  go  to  the  next  step. 

3.  Visit  all  child  nodes  w,  of  v,  and  obtain  probabilities  for  the  sub-functions  fUj  represented  by 

4.  Multiply  the  obtained  probabilities  for  a  sub-function  Ps  (/„.  =  s)  by  the  probability  that  the  component  Xj 
selects  the  sub-function  Pc(x;  =  c). 

5.  Each  function  value  fUj  =  s  at  each  child  node  Uj  becomes  a  function  value  fv  =  s  +  et  at  the  node  v  because 
of  its  edge  value  e,-.  Thus,  the  probabilities  Ps(fUj  =  s)  x  Pc{xi  =  c)  obtained  by  the  step  4  are  added  to 
Ps(fv  =  s  +  ei),  and  they  are  summed  up  (merged)  in  each  function  value  at  v. 

6.  Return  the  merged  probabilities  to  a  parent  node. 

Figure  4.  Proposed  analysis  algorithm  using  EVMDDs. 


NTP(v 3)  =NTP{v  1)  x  1/3  in  a  breadth  first  order.  Sim¬ 
ilarly,  NTP(v4)  =  NTP(vi)/3  +  NTP(v2)/3  +  NTP(v3)/3, 
NTP(v5)  =  2NTP(v2)/3  +  NTP(v3)/3,  and  finally  we 
have  the  node  traversing  probabilities  of  terminal  nodes: 
NTP(0)  =  25/81,  NTP(  1)  =  10/27,  NTP( 2)  =  10/81, 
NTP(3)  =  1/9,  NTP{4)  =  2/27,  and  NTP(5)  =  1/81. 

(End  of  Example) 

B.  Analysis  Method  Using  EVMDDs 

In  an  EVMDD,  a  function  value  is  represented  by  a  sum 
of  edge  values,  rather  than  a  terminal  node.  Thus,  we  cannot 
solve  Problem  1  using  only  node  traversing  probabilities,  and 
another  analysis  method  is  needed. 

Fig.  4  shows  the  proposed  analysis  algorithm.  This  al¬ 
gorithm  visits  each  node  only  once  in  depth  first  order 
starting  from  the  root  node,  and  analyzes  a  sub-function 
represented  by  each  node  recursively.  Probabilities  for  a 
function  represented  by  a  node  can  be  computed  by  merging 
probabilities  for  sub-functions  represented  by  its  child  nodes. 
Thus,  the  algorithm  shown  in  Fig.  4  can  compute  the 
probability  of  each  state  of  a  multi-state  system  correctly 
and  efficiently.  Since  the  algorithm  visits  each  node  only 
once,  its  time  complexity  is  0(Ne),  where  Ne  is  the  number 
of  nodes  in  an  EVMDD. 

Example  4:  Let  us  compute  the  probability  of  each  state 
of  the  multi-state  system  using  the  EVMDD  in  Fig.  3.  As 
with  the  previous  example,  we  assume  that  all  states  of 
each  component  appear  in  the  same  probability  1/3.  First, 
we  have  Ps(fr  =  0)  =  1  at  the  terminal  node  T.  Then,  we 
compute  probabilities  for  a  sub-function  at  the  node  vi. 
Since  this  node  has  two  edges  pointing  to  T  whose  values 
are  1,  we  have 

Ps(fr  =  0)  x  Pc(x4  =  1)  =  1/3, 

Ps(fr  =  0)  x  Pc(x 4  =  2)  =  1/3, and  thus, 


Ps(fVl  =  0+  1)  =  Ps{fr  =  0)  x  Pc(x 4  =  1) 

+Ps(fr  =  0)  x  Pc(x 4  =  2)  =  2/3 
Similarly,  at  the  node  V2,  we  have 

Ps(fv2  =  o  +  o)  =  Ps(fT  =  0)  x  Pc(x 4  =  0)  =  1/3, 

Ps(fv2  =  0  +  1)  =  Psifr  =  o)  x  Pc(x 4  =  1)  =  1/3, and 
Ps(fv2  =0  +  2)=  Ps(fT  =  0)  x  Pc(x 4  =  2)  =  1/3. 

At  V3,  the  probabilities  at  the  terminal  node,  vi,  and  V2 
are  multiplied  by  1/3,  and  are  summed  up.  Thus,  we  have 
Ps(fv 3  =  0)  =  5/9,  Ps(fv 3  =  1)  =  1/3,  and  Ps(fV3  =  2)  =  1/9. 
By  performing  the  same  computation  at  each  node  in  the 
depth  first  order,  we  have  the  following  at  the  root  node: 
P*(f  =  0)  =  25/81,  Ps(f  =  1)  =  10/27,  Ps(f  =  2)  =  10/81, 
Ps{f  =  3)  =  1/9,  Ps(f  =  4)  =  2/27,  and  Ps(f  =  5)  =  1/81. 


Table  I 

MDDS  AND  EVMDDS  FOR  m-STATE  SYSTEMS  WITH  n  3-STATE  COMPONENTS. 


n 

m 

Number  of  nodes 

Computation  time  (//sec.) 

MDD 

EVMDD 

Ratio 

MDD 

EVMDD 

Ratio 

5 

3 

12 

10 

83% 

0.30 

1.20 

406% 

5 

10 

36 

18 

50% 

1.10 

2.52 

230% 

10 

3 

17 

15 

88% 

0.52 

1.85 

355% 

10 

10 

77 

57 

74% 

2.83 

7.84 

277% 

10 

100 

599 

265 

44% 

26.41 

49.27 

187% 

10 

1,000 

4.201 

907 

22% 

231.10 

317.79 

138% 

15 

3 

32 

30 

94% 

1.16 

3.72 

320% 

15 

10 

120 

105 

88% 

4.94 

14.07 

285% 

15 

100 

1,098 

708 

64% 

59.43 

110.72 

186% 

15 

1,000 

9,010 

3,362 

37% 

589.10 

744.59 

126% 

15 

10,000 

70,140 

11,474 

16% 

4,705.00 

4,701.00 

100% 

15 

100,000 

495,224 

62,759 

13% 

65,303.00 

60,901.00 

93% 

n:  Number  of  3-state  components.  m:  Number  of  states  for  systems. 

Ratio:  EVMDD  /  MDD  x  100  (%) 

The  computation  time  is  an  average  time  obtained  by  running  the  same  computation 
a  million  times,  and  dividing  its  total  time  by  a  million. 


Note  that  these  are  consistent  with  the  results  obtained  by 
MDDs  in  Example  3.  (End  of  Example) 

When  a  structure  function  is  monotone  increasing,  the 
number  of  nodes  in  an  EVMDD  Ne  is  smaller  than  for 
non-monotone  increasing  functions,  and  thus,  computation 
time  is  shorter.  Of  course,  the  proposed  method  can  be 
applied  to  nonmonotonic  structure  functions  used  in  some 
applications  [14],  [16]  as  well. 

V.  Experimental  Results 

To  show  the  effectiveness  of  the  proposed  method,  we 
used  various  structure  functions.  Unfortunately,  however, 
benchmark  structure  functions  of  multi-state  systems  were 
unavailable.  Since  structure  functions  are  usually  monotone 
increasing,  we  randomly  generated  Ml -monotone  increas¬ 
ing  functions,  and  used  them  as  structure  functions  for 
experiments  in  this  paper.  The  analysis  algorithms  based  on 
MDDs  and  EVMDDs  are  implemented  using  the  follow¬ 
ing  computer  environment:  CPU:  Intel  Core2  Quad  Q6600 
2.4GHz,  memory:  4GB,  OS:  CentOS  5.7,  and  C-compiler: 
gcc  -02  (version  4.1.2).  Table  I  shows  the  experimental 
results  for  randomly  generated  m-state  systems  with  n  3- 
state  components. 

From  this  table,  we  can  see  that  EVMDDs  have  fewer 
nodes  than  MDDs  for  all  the  functions.  Especially,  as  the 
number  of  states  m  becomes  larger,  EVMDDs  are  much 
smaller  than  MDDs.  We  expect  that  systems  will  become 
more  complex  in  the  future,  and  that  m  would  become  larger. 
Thus,  EVMDDs  whose  size  is  independent  of  the  number 
of  states  are  more  promising.  However,  when  m  is  very 
small,  MDDs  are  better,  since  they  are  small  enough  and 
computation  time  is  shorter.  In  Table  I,  when  m  =  3,  only 
two  terminal  nodes  are  reduced  in  EVMDDs.  Thus,  using 
EVMDDs  for  such  systems  is  not  effective. 

As  for  computation  time,  the  proposed  method  using 
EVMDDs  is  comparable  to  methods  using  MDDs.  There¬ 


fore,  we  can  say  that  EVMDDs  are  suitable  for  compact 
representation  and  efficient  analysis  of  many-state  systems. 

VI.  Conclusion  and  Comments 

This  paper  proposes  an  efficient  analysis  method  of 
multi-state  systems  using  EVMDDs.  The  proposed  analysis 
method  is  somewhat  more  complicated  than  existing  meth¬ 
ods  using  ordinary  MDDs  because  a  state  of  the  system 
is  represented  by  a  sum  of  edge  values.  However,  actual 
computation  time  of  the  proposed  method  is  comparable  to 
methods  using  MDDs  since  the  time  complexity  is  asymp¬ 
totically  proportional  to  the  number  of  nodes  in  an  EVMDD, 
and  EVMDDs  have  fewer  nodes  than  MDDs.  Especially,  for 
systems  with  many  states,  the  proposed  method  is  effective 
because  EVMDDs  are  much  smaller  than  MDDs.  Even  if 
structure  functions  are  not  monotone,  EVMDDs  are  not 
larger  than  MDDs.  Thus,  the  proposed  method  is  effective 
for  a  wide  range  of  structure  functions. 

In  this  paper,  we  used  randomly  generated  Ml -monotone 
increasing  functions  for  our  experiments,  since  benchmarks 
of  multi-state  systems  were  unavailable.  However,  there 
could  be  functions  more  suitable  for  multi-state  systems. 
Thus,  we  will  study  such  functions.  We  will  also  study 
how  to  generate  EVMDDs  directly  from  multi-state  systems 
without  using  MDDs  for  structure  functions. 
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Appendix 

Proof  for  Theorem  1:  Suppose  that  an  MDD  is  parti¬ 
tioned  into  two  parts:  the  upper  and  the  lower  parts  as  shown 
in  Fig.  A.l.  In  this  case,  the  lower  part  represents  /-variable 
multi-valued  functions,  and  the  upper  part  represents  the 
function  that  chooses  one  from  them.  The  upper  part  has  the 
maximum  number  of  nodes  when  it  forms  a  complete  multi¬ 
valued  tree.  The  number  of  nodes  in  the  complete  multi¬ 
valued  tree  is  1  +  r  +  r2  +  . . .  +  Thus,  the  maximum 

number  of  nodes  in  the  upper  part  is 


The  lower  part  has  the  maximum  number  of  nodes  when 
it  represents  all  /-variable  multi-valued  functions.  Since  the 
number  of  all  /-variable  multi-valued  functions  is 

m'J ,  (A. 2) 

that  is  the  maximum  number  of  nodes  in  the  lower  part 
including  terminal  nodes.  From  (A.l)  and  (A. 2),  the  number 
of  nodes  in  the  MDD  is  at  most 

r n~l  -  1  ,./ 

r —  1 

The  number  of  multi-valued  functions  which  can  be  rep¬ 
resented  in  the  lower  part  does  not  exceed  the  number  of 
functions  which  can  be  chosen  by  the  upper  part:  r"  1 . 
Therefore,  we  have  the  relation: 


